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Abstract 

We report on precise hyperspherical-basis calculations of r)NN and rjNNN 
quasibound states, using energy dependent r]N interaction potentials derived 
from coupled-channel models of the Su A* (1535) nucleon resonance. The 
r)N attraction generated in these models is too weak to generate a two-body 
bound state. No r/NN bound-state solution was found in our calculations 
in models where Re a V N < 1 fin, with a V N the i]N scattering length, cover¬ 
ing thereby the majority of A* (1535) resonance models. A near-threshold 
77 AAA bound-state solution, with 7] separation energy of less than 1 MeV 
and width of about 15 MeV, was obtained in the 2005 Green-Wycech model 
where Re a V N ~ 1 fm. The role of handling self consistently the subthreshold 
r]N interaction is carefully studied. 

Keywords: few-body systems, mesic nuclei, forces in hadronic systems and 
effective interactions 


1. Introduction 

The 7]N interaction has been studied extensively in photon- and hadron- 
induced production experiments on free and quasi-free nucleons, and on nu¬ 
clei [l], These experiments suggest that the near-threshold ?/A interaction is 
attractive, but are unable to quantify this statement in any precise manner. 
Nevertheless, r] production data on nuclei provide some useful hints on pos¬ 
sible 7] quasibound states for very light species where, according to Krusche 
and Wilkin (KW) “the most straightforward (but not unique) interpretation 
of the data is that the rjd system is unbound, the r/ 4 He is bound, but that 
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the ? 7 3 He case is ambiguous” [lj. Indeed, the prevailing theoretical consen¬ 
sus since the beginning of the 2000s, based on r)NN Faddeev calculations, 
is that r)d quasibound or resonance states are ruled out for acceptable r]N 
interaction strengths [2j, 0. Instead, the rjd system may admit virtual states 
0, S 0. Searching for reliable few-body calculations of the A = 3,4 rj- 
nuclear systems, we are aware of none for rjNNNN and of only one r]NNN 
Faddeev-Yakubovsky calculation 0], although not sufficiently realistic, that 
finds no ? 7 3 H quasibound state. Rigorous few-body calculations substantiat¬ 
ing the KW conjecture quoted above are therefore called for. The present 
work fills some of this gap, reporting precise calculations of r]NN and of 
7jNNN few-body systems using the hyperspherical basis methodology 0, 
similarly to the calculations reported in Ref. [9i] for the KNN and KNNN 
systems. Particular attention is given in the present work to the subthresh¬ 
old energy dependence of the r/N interaction in a way not explored before in 
rj few-body calculations. 
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Figure 1: Real (left panel) and imaginary (right panel) parts of the r]N cm scattering 
amplitude F V N(y/s) as a function of the total cm energy y/s in five meson-baryon coupled- 
channel interaction models, in decreasing order of Re a V N ■ Dot-dashed curves: GW fllj |; 
solid: CS [l2j; dotted: KSW [h|; long-dashed: M2 [3]; short-dashed: IOV [3|- The thin 
vertical line denotes the rjN threshold. Figure adapted from (l6j - 


Theoretically, the r/N interaction has been studied in coupled-channel 
models that seek to fit or, furthermore, generate dynamically the prominent 
IV* (1535) resonance which peaks about 50 MeV above the r/N threshold. 
Such models result in a wide range of values for the real part of the rjN 
scattering length a V N, from 0.2 fm [l(J to almost 1.0 fm ll|. Most of these 
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analyses constrain the imaginary part Im a V N within a considerably narrower 
range of values, from 0.2 to 0.3 fm. This is demonstrated in Fig. |T] where the 
real and imaginary parts of the qN center-of-mass (cm) scattering amplitude 
F v n(\/s ) are plotted as a function of the cm energy y/s for several coupled 
channel models. The qN threshold, where i^jv(v^th) = is denoted by 
a thin vertical line. We note that both real and imaginary parts of F, qN ( y/s) 
below threshold decrease monotonically in all of these models upon going 
deeper into the subthreshold region, displaying however substantial model 
dependence. This will become important for the q few-body calculations 
reported here. 

Beginning with the pioneering work by Haider and Liu 0. and using 
input values of a v ^ within these specified ranges, several //-nucleus optical- 
model bound-state calculations concluded that n mesons are likely to bind in 

0. In 


sufficiently heavy nuclei, certainly in i3 C and beyond [18|, 119 


21 


the few-body calculations reported here we find no qd quasibound states for 
values of Re a V N as large as about 1 fm. We do find, however, a very weakly 
bound and broad ?/ 3 H-?/ 3 He isodoublet pair for Re a V N ~ 1 fm by solving 
the qNNN four-body problem. 

The paper is organized as follows. In section [2] we construct local energy- 
dependent single-channel potentials v v n that reproduce two of the s-wave 
scattering amplitudes F. qN ( y/s) shown in Fig. Q] GW [U] and CS 12]. In 
section [3] we sketch the hyperspherical-basis formulation and solution of the 
qNN and qNNN Schroedinger equations below threshold using these de¬ 
rived qN potentials and realistic energy-independent NN potentials. Be¬ 
cause of the substantial energy dependence of v v n in the subthreshold re¬ 
gion, a self consistency requirement {sj is applied so that the input energy 
argument of the two-body potential v v n for convergent few-body solutions is 
consistently related to some energy expectation values in the resulting qua¬ 
sibound state. Results are presented and discussed in section [4] followed by 
a brief summary and outlook section O 


2. Construction of r]N effective potentials 

We seek to construct energy-dependent local qN potentials v v n that re¬ 
produce the qN scattering amplitude F v jy(^/s) below threshold in given mod¬ 
els, e.g. from among those shown in Fig. [TJ For convenience, the energy argu¬ 
ment E introduced in this section is defined with respect to the qN threshold, 
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E = y/s — and should not be confused with the binding energy of the 

r]NN and rjNNN few-body states studied in subsequent sections. 

We define v v n in the form 

v V N(E;r) — —^—b(E)p A (r), (h = c = 1) (1) 

with /j, v jsf the reduced r/N mass and where p\ is a Gaussian normalized to 1: 

Mr) = (2^) exp ■ (2) 

A is a scale parameter, inversely proportional to the range of v v n. Its physi¬ 
cally admissible values are discussed in subsection 12.21 below. Two represen¬ 
tative values are used here, A=2 and 4 fm -1 . For a given value of A, one 
needs to determine the energy-dependent strength parameter b(E) of v v n, as 
described in the following subsection 12.11 


2.1. Solution 

Given a specific value of the scale parameter A, the two-body s-wave 
Schroedinger equation 

- 7 T— U '\ r ) + v vN (E] r)u(r ) = Eu(r) (3) 

Zli v N 

is solved for energies above (E > 0) and below (E < 0) threshold. The radial 
wavefunction u{r) satisfies the boundary conditions 


u{r = 0) = 0, u{r —> 00 ) oc r(cos5 0 jo(kr) — sin^ 0 n 0 (kr)), (4) 


where k = yj2 /i^nE, jo and no are spherical Bessel and Neumann functions, 
respectively, and So(E) is the complex s-wave phase shift derived by imposing 
these boundary conditions on the wave-equation solution. Above threshold, 
the wave number k is real and taken positive. Below threshold, k = ik with 
n > 0. The scattering amplitude F is then given by 


F v n(E) — 


k( cot 5 0 — i ) 


(5) 


This procedure was used in Ref. 23] for constructing effective KN potentials 
below threshold. In the present case, the subthreshold values of the complex 
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strength parameter b(E) in Eq. (JT]) were fitted to the complex phase shifts 
6(E) derived from subthreshold scattering amplitudes F vN (E ) in several of 
the coupled-channel models of Fig. [I] This is shown for the GW 0 and 
CS [|l2| models in Fig. [2] using two values of the scale parameter A=2 and 
4 fm 1 for GW and just one value A = 4 fm -1 for CS. The curves b(E) are 
seen to decrease monotonically in going deeper below threshold, except for 
small kinks near threshold that reflect the threshold cusp of Re E v n(E — 0) 
in Fig. [0 Comparing models GW and CS for the same scale parameter 
A = 4 fm -1 , one observes larger values of b(E) in model GW than in CS, for 
both real and imaginary parts below threshold, in line with the larger GW 
subthreshold amplitudes compared with the corresponding CS amplitudes. 
We note furthermore that Im b(E )<CRe b(E) in both models by almost an 
order of magnitude, see Fig. [2l which justifies treating Im v v n perturbatively 
in the applications presented below. 



Figure 2: Real (left panel) and imaginary (right panel) parts of the strength parameter 
b(E) of the r/N effective potential (fljl. for subthreshold energies E < 0, obtained from the 
scattering amplitudes F^ [llf and jli} shown in Fig. |T] Two choices of the scale 
parameter A are made for GW, both resulting in the same F^(E), and just one for CS. 


To demonstrate the extent to which the energy dependence of b(E) is 
essential, we compare in Fig. [3] the GW subthreshold amplitude from Fig. |T] 
which is also generated here using the b(E) potential strength of Fig. [2] for 
A = 4 fm -1 , to the amplitude marked gw which was calculated using a fixed 
threshold value b(E = 0). This latter amplitude is seen to decrease too slowly 
beginning about E ~ —7 MeV. Obviously, an energy-independent single- 
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Figure 3: Real (left panel) and imaginary (right panel) parts of the subthreshold amplitude 
( E) (solid curves marked GW) from Fig. [I] also generated from an energy dependent 
potential ( E ) with A = 4 fm 1 , Eqs. ( 11121 ) . compared to the amplitude (dashed curves 

marked gw) generated from (E = 0). 


channel potential v v n fails to reproduce the subthreshold energy dependence 
of the GW coupled-channel scattering amplitude F^(E). 

2.2. Choice of scale 

It is appropriate at this point to address the model dependence intro¬ 
duced in //-nuclear few-body calculations by the choice of the scale parameter 
A made in constructing v v n, Eqs. (jTj) and ([2]) . A is often identified with the 
momentum cutoff used to renormalize divergent loop integrals in on-shell 
EFT V*(1535) models [14, [15]. In separable-interaction coupled channel 


models, however, the momentum cutoff is replaced by fitted Yamaguchi form 
factors ( q 2 + A 2 ) -1 with a momentum-space range parameter A, the Fourier 
transform of which is a Yukawa potential exp(—Ar)/r with r.m.s radius iden¬ 
tical to that of the Gaussian potential shape (J2])- Values of A from three such 
V*(1535) models, including the two used in the present work [111 . 12], are 
listed in Table [H 

Inspection of Table [j] reveals a broad range of values that A may assume, 
starting with A k 3 fm -1 . The relatively high value in the third column 
is rather exceptional for meson-baryon separable models. Given this broad 
spectrum of values spanned for A, we chose two values A = 2 and A = 4 fm -1 
to study the model dependence of our //-nuclear few-body calculations. The 
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Table 1: The r/N momentum scale parameter A from several _/V*(1535) separable models. 


Ref. 

[10,13J 

[11] 

112 

A (fm -1 ) 

3.9 

3.2 

6.6 


higher value, A = 4 fm -1 , corresponds to a Gaussian exp(— r 1 2 /R 2 ) spatial 
range R = 2/A = 0.5 fm, a value which is very close to R = 0.47 fm taken 
from the systematic EFT approach in Ref. [23j] and used in our A-nuclear 
few-body calculations |9!]. As argued there, choosing smaller values for R, 
namely larger values than 4 fm -1 for A, would be inconsistent with staying 
within a purely hadronic basisQ 

In the Introduction section we loosely identified the strength of the r]N 
interaction with the size of the real part of its threshold scattering amplitude, 
Re UriN < 1 f m - However, in terms of the interaction potentials v v n that enter 
our few-body calculations, a given value of Re a VJ \ does not rule out a broad 
spectrum of spatial ranges, or equivalently momentum scale parameters A, 
as demonstrated in Fig. [2] A model dependence is thereby introduced into 
our few-body calculations, summarized by stating that the larger the r)N 
scale parameter A is, the larger is the rj separation energy, provided it is 
quasibound. This lack of scale invariance hints towards the necessity of 
including three-body forces, as is expected from an EFT point of view [25|. 
Such three-body forces amount to adding a new free parameter determined 
by tuning it to some r) few-body experimental data. 


1 The effective energy-dependent KN potential Vf <N constructed by Hyodo and Weise 

[ 23 I] reproduces the KN — 7 tE coupled-channel scattering amplitude which is the one es¬ 
sential for generating dynamically the A*(1405) resonance. In that case, the choice of A 
must ensure that the K*N channel that couples strongly to KN via normal pion exchange 
is kept outside of the model space in which v^ N is valid. This argument leads to a choice 
of A = p m in(KN — > K*N ) = 552 MeV/c or 2.8 fm -1 , corresponding to a Gaussian spatial 
range of R = 0.71 fm. In a somewhat similar reasoning Garzon and Oset [24[ recently 
argued for extending the EFT description of the N* (1535) resonance to include the pN 
channel which couples strongly to the already included nN channel, although not to r/N. 
Identifying A with the minimum momentum needed to excite the irN system to pN, we 
obtain A = p m i n (nN —t pN) = 586 MeV/c or 3.0 fm -1 . 
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3. 77 -nuc l ear hyperspherical-basis formulation and solution 

The hyperspherical-basis formulation of meson-nuclear few-body calcu¬ 
lations was initiated in Ref. {§] for K mesons. Here we sketch briefly the 
necessary transformation from K mesons to 7] mesons. The iV-body wave- 
function (N = 3,4) in our case consists of a sum over products of isospin, 
spin and spatial components, antisymmetrized with respect to nucleons. In 
the spatial sector translationally invariant basis functions are constructed in 
terms of one hyper-radial coordinate p and a set of 3 N — 4 angular coordi¬ 
nates [fl at], substituting for N — 1 Jacobi vectors. The spatial basis functions 
are of the form 

*»,i*i(p,n K ) = J^(p)3' [ 5](n N ), (6) 

where (p) are hyper-radial basis functions expressible in terms of La- 
guerre polynomials and are hyperspherical-harmonics (HH) func¬ 

tions in the angular coordinates expressible in terms of spherical har¬ 
monics and Jacobi polynomials. Here, the symbol [K\ stands for a set of 
angular-momentum quantum numbers, including those of L 2 , L z and K 2 , 
where K is the total grand angular momentum which reduces to the total 
orbital angular momentum for N = 2. The HH functions Tyy are eigen¬ 
functions of K 2 with eigenvalues Ii(K + 3 N — 5), and p K y^ K i are harmonic 
polynomials of degree K [§j. 

For the NN interaction we used two forms, the (Minnesota) MN central 
potential [26j] and the Argonne AV4’ potential [27j derived from the full 
AV18 potential by suppressing the spin-orbit and tensor interactions and 
readjusting the central spin and isospin dependent interactions. In s-shcll 
nuclei the AV4’ potential provides an excellent approximation to AV18 which 
pseudoscalar mesons, such as the r) meson, are unlikely to spoil, recalling that 
their nuclear interactions cannot induce S yy D mixing beyond that already 
accounted for by the NN interaction!! AV4’ and MN differ mostly in their 
short-range repulsion which is much stronger in AV4’ than in MN. 

For the r]N interaction we used the energy-dependent local potential 
Re v v n introduced in Sect. [21 In order to distinguish the energy E of the 


2 This was demonstrated in K nuclear cluster calculations Q, see the discussion of 
Table 1 therein, where the K(NN)i—o 4.7 MeV binding energy contribution to the full 
15.7 MeV binding energy of (KNN)i-u 9 calculated using AV4’ is short by only 0.2 MeV 
from that in a comparable calculation [28| using AV18. 



few-body system from the energy argument of v vN , the latter is replaced by 
Sy/s = y/s — yfsth from now on. Following Eq. (5) in [oj], the subthreshold 
energy argument 5^/s of v v n, is chosen to agree self-consistently with 

= -§ - - £„ (^) W- (?) 


where £ at ( ?7 ) = rnN( v )/(mN + m v ), T v is the rj kinetic energy operator in the 
total cm frame, T^-.n is the pairwise NN kinetic energy operator in the NN 
pair cm frame, B is the total binding energy of the //-nuclear few-body system 
and B v is the rj “binding energy”, B v = —E^ = —— F/at)|'F), where 
Hn is the Hamiltonian of the purely nuclear part in its own cm frame and the 
total Hamiltonian H is evaluated in the overall cm frame. I 11 the limit A 1, 


Eq. ((7|) agrees with the nuclear-matter expression given in Refs. 2l|, 22] for 
use in calculating -//-nuclear quasibound states. It provides a self-consistency 
cycle in 77 -nuclear few-body calculations by requiring that the expectation 
value (5i/s) derived from the solution of the Schroedinger equation agrees 
with the input value 5^/s used in v v m- Since each one of the four terms on the 
r.h.s. of (JTJ) is negative, the self consistent energy shift 5y/s s , c , is necessarily 
negative, with size exceeding a minimum nonzero value obtained from the 
first two terms in the limit of vanishing 7 / binding. 

The potential and kinetic energy matrix elements for a given rj -nuclear 
state with global quantum numbers /, L, S , J 71 " were evaluated in the HH ba¬ 
sis. The NN and 7/ IV interactions specified above conserve / = In, S = Sn 
and L. Since no L ^ 0 7/-nuclear states are likely to come out particle stable, 
our calculations arc limited to L = 0. The deuteron in this approximation is 
a purely 3 S± state. Suppressing Imr^, the g.s. energy E g s . was calculated 
in a model space spanned by HH basis functions with eigenvalues K < /i max . 
Self-consistent calculations were done for a/s ranging from the r]N thresh¬ 
old down to 30 MeV below. Self consistency in 5^fs was reached after a 
few cycles. Good convergence was achieved for values of /l max « 20 — 40. 
Asymptotic values of E g s . were found by fitting the constants C and a of the 
parametrization 


Ei^I\. max) fr'g.s. T C exp( OiE max) 


( 8 ) 


to values of E( K max ) calculated for sufficiently high values of K max . The 
accuracy reached is better than 0.1 MeV in both the three-body and the 
four-body calculations reported here. 
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The conversion width T was then evaluated through the expression 

T = -2(T g . s . |ImV^v|tf g . B .) , (9) 

where V v n sums over all pairwise rjN interactions. Since |Im V^jvl |Re V v n\, 
this is a reasonable approximation for the width. 


4. Results and discussion 

Results of r)NN and r/NNN bound-state hyperspherical-based calcula¬ 
tions for the GW rjN interaction, with Re a v n almost 1 fm, are given in 
this section. The weaker CS rjN interaction is found too weak to generate 
bound-state solutions. 


J h l. r]NN calculations 

No 1 = 0, = 1 “ ijd bound state solution was found for the rjNN three- 

body system using the MN NN potential j26[ and the GW nj rjN effective 
potential with a fixed strength b(5^/s = 0), see Fig. El for either choice A = 2 
or 4 fna 1 of the scale parameter under study. It was found that b(5y/s = 0) 
in the GW model needs to be multiplied by 1.1 for A = 4 fm -1 and by 1.3 for 
A = 2 fna 1 in order to generate a 1“ r/NN weakly bound state, with overall 
binding energy of —2.219 and —2.264 MeV, respectively, within three-body 
calculations that use a fixed rjN interaction strength b(5y/s = 0). Recall 
that the MN deuteron binding energy is Ed = —2.202 MeV. There is no rjd 
bound state also in the r]N CS jlij model, judging by the CS/GW relative 
strengths of b(5y/s). 

Given that the rjN interaction is too weak to bind the 1 = 0, J n = 1“ 
r/NN state in which the 3 S'i NN (deuteron) core configuration is bound, 
the unbound 1 S'o NN core configuration in the 1 = 1, J n = 0~ rjNN state 
certainly cannot support a three-body bound state. This holds so long as 
the 1“ state is unbound and also for a certain range of larger rjN potential 
strengths that make the 1“ bound. This situation is reminiscent of the 
ANN system which is known to have one 1 = 0 bound state in which the A 
hyperon is bound to a deuteron core, but no I = 1 ANN bound state, see 


e.g. Ref. 29 


Our negative results rule out any rjd bound state, practically in all dynam¬ 
ical models of the V*(1535) resonance where the rjN interaction is coupled 
in, and are consistent with similar conclusions reached in Refs. [2, [3, 4J, |5|, |6|. 
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This holds also upon replacing the MN NN interaction 26] by the AV4’ NN 
interaction [27| in our r]NN calculations. In fact, somewhat larger rjN in¬ 
teraction multiplicative factors are then required to reach the onset of r]NN 
binding compared to those specified above. Applying the self-consistency 
requirement discussed in Sect. [3] to the r t]NN calculation, and recalling the 
decreased strength b(5y/s) in the rjN subthreshold region, see Fig. [21 would 
only aggravate the failure to generate a three-body rjNN bound state. 


4-2. r]NNN calculations 

Four-body rjNNN calculations were made using the MN [26| and the 
AV4’ [i?} NN potentials, and the GW 11] and CS 12] energy-dependent rjN 


potentials from Sect. [21 Based on r/ 3 H and rf He, and with the leading 3N con¬ 
figuration given by In = Sn = \ and I7 N = 0 + , the quantum numbers of the 
calculated rjNNN state are I = S = L = 0 and J n — \ . The 3N bind¬ 
ing energy (disregarding the Coulomb interaction in the case of 3 He) within 
our hyperspherical-basis calculation is —8.38 MeV for MN and —8.99 MeV 
for AV4\ Starting with the rjN GW model, with Re a V N = 0.96 fm, and 
using the corresponding v v n from Sect. [ 2 ] with energy independent threshold 
strength b(5y/s = 0) for A = 4 fm -1 , a four-body r/NNN bound state was 
found with ij separation energy E s ’ a between 2 to 3 MeV, as listed in Ta¬ 
ble [21 We then applied a self consistency procedure by doing calculations 
with several given values of strength b(5y/s), requiring that the expectation 
value (5y/s) evaluated by Eq. (J7J) from the obtained solution agrees with the 
input value of the subthreshold energy Sy/s argument of the strength b(5y/s) 
used in the calculation. This resulted in considerably reduced values of less 
than 1 MeV for the rj separation energy E^ c sep which are listed in Table El 
together with the corresponding r]NNN binding energies E^\ Also listed in 
the table are the self consistent values 5^/s s . c . and the self-consistency reduc¬ 
tion factors x s . c . = b(5y/s s , c ,)/b(5y/s = 0). No rjNNN bound-state solutions 
were found using vffl self consistently for A = 2 fm -1 . 


Table 2: Results of yNNN quasibound-state self-consistent calculations using the r)N 
model GW [Tl]. Energies and widths are given in MeV. 


NN int. 

E(NNN) 

ZT’no s.c. 

Kb- 

by/ $s.c. 

^S.C. 

77 ’S.C. 

Kb. 

771 s.c. 

^r) sep. 

■p S.C. 

g.s. 

MN 

-8.38 

-11.26 

-13.52 

0.914 

-9.33 

0.95 

13.52 

AV4’ 

-8.99 

-11.33 

-15.83 

0.895 

-9.03 

0.04 

15.75 
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Figure 4: The rjNNN g.s. energy E g . s . (solid curves) and the expectation value (Sy/s) 
(dashed curves) from Eq. iff]), calculated using the NN potentials MN (red) and AV4’ 
(blue), are shown as a function of the energy argument Sy/s used for the input tVy V . The 
dashed horizontal line marks the NNN ( 3 H) g.s. energy —8.48 MeV and the dashed 
diagonal line marks potentially self consistent solutions satisfying (Sy/s) = Sy/s. The 
dashed vertical lines mark the intersection of the dashed diagonal line with the (Sy/s) 
dashed curves, thereby fixing the self-consistent values Sy/s s . c .- 


In order to demonstrate how the self consistency procedure works we 
plotted in Fig. 0 the 7]NNN g.s. energy E gs ^ and expectation value (Sy/s), 
calculated as a function of the subthreshold energy Sy/s argument of the input 
r/N potential v rjN in both NN potential models. The difference between the 
.Eg.s. curves, using MN or AV4’, is a fraction of MeV for any given input 
value Sy/s and is hardly noticeable in the figure. The difference between the 
corresponding (Sy/s) curves amounts to a few MeV at each value of Sy/s and 
is clearly visible in the figure, leading to self-consistency values 5^/s sx . which 
differ from each other by more than 2 MeV (marked by the dashed vertical 
lines). The corresponding self consistent values of E gmSm are much closer to 
each other (marked by the thin dashed horizontal lines). The self consistency 
procedure is applied in the figure by looking for the intersection of the dashed 
diagonal line, locus of (5\/i) = Sy/s, with each of the (Sy/s) dashed curves. 
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Applying a similar self-consistency procedure to the weaker CS ijN inter¬ 
action, rather than to the GW r]N interaction used above, no ijNNN bound 
state solution was found. With AV4’ for the NN interaction, this holds even 
upon using the threshold energy value in v^. With the MN NN interaction 
and for the choice A = 4 fm _1 , a bound-state solution is found for small val¬ 
ues of the input energy 5^/s, disappearing at — 5y/s ~ 2.8 MeV which is way 
below the minimum value of —6y/s required in the limit of E v sep —y 0. We 
conclude that the CS r]N interaction is too weak to provide self consistently 
r]NNN bound states. 

Finally, the r]NNN width ~ 15 MeV listed in the last column of 
Table [2] was calculated using Im b(Sy/s s , c ,) in forming the integrand Im V v n in 
Eq. ([9]) . This width is about three times larger than the widths evaluated self 
consistently using optical model methods across the periodic table within the 
rjN GW model [21], Some explanation of this difference is offered noting that 
the magnitude of the downward energy shifts 5y/s s , c . effective in those works 
is considerably larger by factors of two to three than the f»15 MeV found in 
the present r/NNN calculations, reflecting the denser nuclear environment 
encountered by the /? meson as it becomes progressively more bound in the 
calculations of Ref. [21], Recalling the steady decrease of the r]N absorptivity 
Im F^n in Fig. |Tj upon moving deeper into subthreshold energies, a factor 
of two to three difference could be anticipated in favor of relatively small 7] 
widths in heavier nuclei. 

5. Summary and outlook 

Precise hyperspherical-based few-body calculations were reported in this 
work to explore computationally whether or not t] mesons bind in light nu¬ 
clei. To this end, complex energy-dependent local effective tjN potentials 
v v n were constructed, for subthreshold energies relevant to 77 mesic nuclei, 
from coupled channel tjN scattering amplitudes in several models connected 
dynamically to the V*(1535) resonance. The scale dependence arising from 
working with an effective v v n was studied by using two representative values 
for the momentum scale, A = 2,4 fm -1 . Noting that Im v v m Re v v n, 
only the real part of v v n was used in the bound-state calculations, with a 
related error estimated as less than 0.2 MeV, added to an estimated 0.1 MeV 
calculational error. The width of the bound state, making it into a quasi¬ 
bound state, was deduced from the expectation value of Im v v n summed on 
all nucleons. 
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No rjNN quasibound states were found for any of the two scale parameters 
chosen in models where the real part of the rjN threshold interaction satisfies 
Re a V N < 1 fm, in agreement with deductions made in several past few-body 

calculations of the rjd scattering length 0 , 00 , 0 , 0 . It is unlikely that the rjd 

system can reach binding upon increasing moderately the momentum scale 
parameter A. 

For rjNNN , essentially the ? 7 3 H and p 3 He isodoublet of r) rnesic nuclei, a 
relatively broad and weakly bound state was found with p separation energy 


of less than 1 MeV using the GW rjN interaction model |ll| where Re a V N is 
almost 1 fm. This holds for the larger of the two values of momentum scale 
parameter, A = 4 fm -1 , studied here, whereas no bound state was obtained 
upon using the smaller value of A = 2 fm -1 . The e nerg y dependence of g gw 


pN 


treated here within a self consistent procedure [21, 22], played an important 
role by reducing the calculated binding energy by about 2 MeV from that 
calculated upon using the rjN threshold energy value in tSV For such halo- 
like //-nuclear quasibound states, the neglect of Im v v n in the bound-state 
calculation requires attention. In the case of the GW rjN effective interaction, 
we estimate the repulsion added by reinstating Im vffl to second order to be 
roughly < 0.2 MeV, eliminating thereby the very weakly bound rjNNN state 
calculated here using the AV4’ NN potential, but not the weakly bound one 
calculated using the MN NN potential. It is worth noting that the only 
other few-body rjNNN study known to us 0 deduced from their calculated 
r/ 3 H scattering length that no quasibound state was likely. However, the 
strength of the rjN interaction tested in these calculations was limited to 
Re a V N = 0.75 fm, short of our upper value of approximately 1 fm. 

In conclusion, recalling the KW conjecture [lj] quoted in the Introduction, 
it is fair to say that the present few-body calculations support the conjecture’s 
first and last items, namely that “the rjd system is unbound” and “that 
the r/ 3 He case is ambiguous”. Accepting that the strength of the two-body 
rjN interaction indeed satisfies Re a V N < 1 fm, which is much too weak 
to bind the rjN system, a persistent theoretical ambiguity connected with 
choosing a physically admissible range of values for the rjN scale parameter 
A is demonstrated by our few-body calculational results, particularly for the 
four-body rjNNN system. By choosing a considerably larger value of A than 
done here one could bind solidly this system. To remove this ambiguity, 
many-body repulsive interactions involving the rj meson need to be derived 
and incorporated within few-body calculations. 

In future work we hope to extend our r/NNN calculations also by applying 
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methods of complex scaling that should enable one to follow trajectories 
of 5-matrix quasibound-state poles and look also for other types of poles 
such as virtual-state poles or resonance poles, all of which affect to some 
degree the threshold production features of rj mesons in association with 
3 He. Furthermore we hope to initiate a precise and realistic calculation of 
the r r/NNNN system in order to test the middle item in the KW conjecture, 
namely that “? 7 4 5 6 7 8 9 10 He is bound”. 

Acknowledgments 

We thank Jiff Mares for carefully reading this manuscript and making 
useful remarks. This work was supported in part (NB) by the Israel Science 
Foundation grant 954/09, and in part (EF, AG) by the EU initiative FP7, 
Hadron-Physics3, under the SPHERE and LEANNIS cooperation programs. 

References 

[1] B. Krusche, C. Wilkin, Prog. Part. Nucl. Phys. 80 (2015) 43. 

[2] A. Deloff, Phys. Rev. C 61 (2000) 024004. 

[3] H. Garcilazo, M.T. Pena, Phys. Rev. C 61 (2000) 064010, 63 (2001) 
021001 (R). 

[4] A. Fix, H. Arenhovel, Eur. Phys. J. A 9 (2000) 119, Nucl. Phys. A 697 
(2002) 277. 

[5] S. Wycech, A.M. Green, Phys. Rev. C 64 (2001) 045206 

[6] H. Garcilazo, Phys. Rev. C 67 (2003) 067001. 

[7] A. Fix, H. Arenhovel, Phys. Rev. C 66 (2002) 024002. 

[8] For a recent review of the hyperspherical basis methodology and other 
ones in few-body calculations, see W. Leidemann and G. Orlandini, 
Prog. Part. Nucl. Phys. 68 (2013) 158. 

[9] N. Barnea, A. Gal, E.Z. Liverts, Phys. Lett. B 712 (2012) 132. 

[10] N. Kaiser, T. Waas, W. Weise, Nucl. Phys. A 612 (1997) 297. 


15 



[11] A.M. Green, S. Wycech, Phys. Rev. C 71 (2005) 014001, 72 (2005) 
029902(E). 

[12] A. Cieply, J. Smejkal, Nucl. Phys. A 919 (2013) 46. 

[13] N. Kaiser, P.B. Siegel, W. Weise, Phys. Lett. B 362 (1995) 23. 

[14] M. Mai, P.C. Bruns, U.-G. Meifiner, Phys. Rev. D 86 (2012) 094033. 

[15] T. Inoue, E. Oset, M.J. Vicente Vacas, Phys. Rev. C 65 (2002) 035204. 

[16] A. Gal, E. Friedman, N. Barnea, A. Cieply, J. Mares, D. Gazda, Acta 
Phys. Pol. B 45 (2014) 673. 

[17] Q. Haider, L.C. Lin, Phys. Lett. B 172 (1986) 257; L.C. Liu, Q. Haider, 
Phys. Rev. C 34 (1986) 1845. 

[18] Q. Haider, L.C. Liu, Phys. Rev. C 66 (2002) 045208. 

[19] C. Garcia-Recio, T. Inoue, J. Nieves, E. Oset, Phys. Lett. B 550 (2002) 
47. 

[20] D. Jido, H. Nagahiro, S. Hirenzaki, Phys. Rev. C 66 (2002) 045202. 

[21] E. Friedman, A. Gal, J. Mares, Phys. Lett. B 725 (2013) 334. 

[22] A. Cieply, E. Friedman, A. Gal, J. Mares, Nucl. Phys. A 925 (2014) 126. 

[23] T. Hyodo, W. Weise, Phys. Rev. C 77 (2008) 035204. 

[24] E.J. Garzon, E. Oset, Phys. Rev. C 91 (2015) 025201. 

[25] P.F. Bedaque, U. van Kolck, Annu. Rev. Nucl. Part. Sci. 52 (2002) 339. 

[26] D.R. Thompson, M. LeMere, Y.C. Tang, Nucl. Phys. A 286 (1977) 53. 

[27] R.B. Wiringa, S.C. Pieper, Phys. Rev. Lett. 89 (2002) 182501. 

[28] A. Dote, T. Hyodo, W. Weise, Nucl. Phys. A 804 (2008) 197, Phys. Rev. 
C 79 (2009) 014003. 

[29] A. Gal, H. Garcilazo, Phys. Lett. B 736 (2014) 93. 


16 



